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Abstract 

The discrete Schrodinger equation with potential belonging to F2 is solved ex- 
plicitly. On this base the associated (1 + l)-dimensional cellular automaton is 
examined and corresponding set of integrals of motions is constructed. 

1 Introduction 

Cellular automata are dynamical systems in the discrete space-time with values belong- 
ing to some finite field. For example, in the simplest situation the dependent variable q l n 
takes values and 1, i.e., belongs to F 2 for any n, t E Z, where n and t are considered 
as space and time variables correspondingly. Initial data are given as a set of zeroes and 
ones numbered by n £ Z and time evolution is a law that enables us to construct 
also as a set of zeroes and ones if q\ L is known. Cellular automata attract essential inter- 
est in the literature because of the wide range of their applications in different sciences, 
from physics to biology, from chemistry to social sciences. Detailed references for these 
and many other applications can be found in [1-5]. These automata are also subject to 
intensive mathematical study, see for example [6-18], but nevertheless the mathematical 
theory of cellular automata is far from completeness. In some cases (see [4]) these dy- 
namical system demonstrate stochastic behavior, in other cases they behave as solitonic 
systems [7-18]. The situation becomes even more intriqueing if to take into account that, 
at least, two families of the cellular automata have Lax pairs. Nevertheless, the problem 
of integrability of such cellular automata is still opened as these Lax representation are 
satisfied only on corresponding finite fields as well. 

Above mentioned families of the cellular automata that have associated Lax repre- 
sentations were considered in [12-15]. All of them are of so called filter type [6]. In 
works [12,15] these were automata with Lax operator L to be the known Ablowitz-Ladik 
system [19]. In work [13] L was chosen to be the discrete (1-dimensional) Schrodinger 
operator, that is also known in the literature [20,21]. As was mentioned in [12,13,15] (and 
also discussed below) known results on either Ablowitz-Ladik, or discrete Schrodinger 
equation cannot be applied to the case where potential belongs to a finite field as in this 



case some of the Jost solutions do not exist. Correspondingly, this unable us to introduce 
the spectral data, so the standard scheme [22-24] of the inverse spectral transform fails 
for this kind of problems requires essential modification. 

In this note we demonstrate that, on the other side, in consideration of the dis- 
crete problems some essential simplifications occur. In particular we show that the Jost 
solutions of the discrete Schrodinger equation (the existing ones!) can be constructed 
explicitly. We apply result of this construction for the simplest example of the cellular 
automata introduced in [13], 

Jn 1 =Ql + C-^n+l + <£t-l<L+2, mod 2. (1.1) 

Here dependent variable q l n £ {0, 1} for any values of two independent variables n,t 6 Z. 
Lax representation for this equation is given in [13] as 

L t+1 A = AL\ mod 2, (1.2) 

where matrices L t mn and A mn are equal to 



L m,nH = S m,n+1 + i 1 + <D 5 m,n-X ~ ( w + S n 



(1.3) 



^m,n — ^m,n ^m+2,nQm-l1 m +2- (1-4) 

It was also shown in [13] that if we consider q l n at some moment of t to be finitely 
supported then and so on also have finite support. More exactly, let us introduce 
support of sequence at moment t, 

IC t = {k t 1 <...<k t Nt } = {kE Z\4 = 1}, (1.5) 

so that k\ and k N t are bottom and upper borders of the support and N < oo. Then 
borders of the support are preserved under time evolution, 

k-j^ — k^, k^ t+1 = kjyt, (1-6) 

and, moreover, support fC is decomposed into "islands" — subsets of support separated 
by three or more consequent sites occupied by zeros. These islands evolve independently, 
in particular, their borders are also integrals of motion. Simple example of time evolution 



of cellular automata (1.1) is presented in Fig. 1, where empty boxes correspond to zero 
sites and black boxes correspond to the sites occupied by units. We presented here the 
case of two islands. It is clear that time evolution is reduced to the motion of single and 
double defects inside islands, and that it is enough to consider the case where K} is a 



single island. Finally, let us notice that equation (1.1) is invariant with respect to the 
substitution q l n — ► qZ^ 1 , so this cellul ar a utomaton is time reversible. 



Spectral problem for the operator (1.3) 



i4,nHV4H = o, (i.7) 



here also has to be understood modulo 2. This means that for any n and t solution 
tpn(w) is formal Laurent series in w 6 C with coefficients belonging to F2. Omitting for 
simplicity subscript t we rewrite (L7) by means of ( [L3| ) as 



V'm-l + (! +9m)VWl = ( W + — ) % l J m, mod 2. 
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This is known discretization of the famous Schrodinger equation and for generic q m 
it can be studied in analogy with the standard Schrodinger equation, see [20,21]. We 
already mentioned that it is just modulo 2 condition that prevents from the use of the 
corresponding results of [20,21] here. Indeed, for all m £ K. (i.e., q m = 0) we can 
determine f/> m +i if V'm and ipm -i are known. On the other side, if m 6 /C (i.e., q m = 1), 
then VVi+i drops out from (1.8) and 7/> m and ip m -i cannot be chosen independently. This 
means that in all cases where K, is not empty equation (1.8) can be solved by means of 
the sweeping from the left, but cannot be swept from the right. 



2 Solution of the discrete Schrodinger 
equation 

In order to construct solutions of the discrete Schrodinger equation let us mention hrst 
that in the case where q m € ¥2 we have 

(l + 5 J(mod 2) = 1 - q m , (2.1) 

so instead of ( |l.8| ) we start with construction of solutions of the infinite system 



(2.2) 



where q m £ {0, 1} for any m £ Z, w € C and condition that (2.2) is satisfied by mod 2 
is omitted. It is clear that such solutions of ( |2.2| ) resolves ( jl.Sj ) in the sense mentioned 
in Introduction. 

If m < ki or m > fcjy + 1 (see (1.5)) we have two explicit solutions of (2.2): w m and 
w~ m . Thus for such m general solution of (2.2) is given as linear combination of these 
two solutions. Let 

x m (w)=w m TP m (w), (2.3) 



that due to (2.2) obeys equation 



(1 + z)x rn + z{l - q m )x m+1 = 0, 



(2.4) 



where for simplicity we denoted 

z = w~ 2 . (2.5) 

In what follows we use the same notations, say x m (z) or x m (w), for functions of cither 
z, or w, related by means of this substitution. It is clear that if ip m (w) solves (2.2) then 

V'mO) =ip m {w~ 1 ) (2.6) 
Let us introduce x m in analogy with ( |2.3| ), then 

x m {z) = z-^x^z- 1 ) (2.7) 



gives another solution of (E 



gives another solution of (|]j). 

Mentioned in Introduction problem of the study of the spectrum (2.4) is reflected in 
the existence of the trivial solution 



b 0,m 



» 



2 z^ N ~ 7n 

1-l/z ' 

0, 



m > k N , 

771 < k N . 



(2.8) 
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This solution is invariant (up to m-independent factor) with respect to transforma- 
tion <$LT\). 

The Jost solution is defined by condition limm^+oo x m — > 1. Then by means of ( |2.4| ) 
we have 

(2.9) 



x m — 1, m > fcjv, 



and 



x m _ x - x m = zq m + o(z), \z\ -> 0, (2.10) 

For m < k ± — 2 this solution is linear combination of 1 and z' m . It can be equivalently 
determined by means of the "integral" equation 

2 gk—m 



El — ^ 
lb I _ I 



-1/z 



(2-11) 



In order to resolve this equation we mention first that due to (1.5) it can be rewritten as 

(2.12) 



^ ykj m - \ _ y 

1 ~ Y-~z Xk * +1 ' 

j:kj>m 



Thus we see that all x m for m — l^K, are determined by those values for which m — 1 G K,. 
Substituting in ( 2.12 ) m by ki + 1, ki € calK. we get for these values the following linear 
algebraic system 



jv 



2; J I z 



E i_ z %+i = 1 > l = l,...,N. 



3=1 



Let us denote 



N 



a i( z ) = 1 + Y^E^+i' 



3=1 



i w 



(2.13) 

(2.14) 
(2.15) 



where Z = 1, . . . , AT. These functions are not independent as because of ( 2.13 ) we have 



a, +z- k '- 1 b l = 0, l=l,...,N, 



so that (2.12) takes the form 



1, m > k N + 1 

a ;+1 (l-z fc <+i-™ +1 ), fc, +1 >m>fc, + l' 



where Z = 0, . . . , N — 1, and we introduced fcp = — oo. 



On the other side, considering difference ct;+i — a; by (2.14) we get 
x kl+1 = (1 - l/z)(a l+1 - a t ), I = 1,...,N- 1. 



From ( 2.14 ) at Z = iV and the upper line of ( 2.17 ) we get that 

1 



'■N 



1 - Z 



(2.16) 



(2.17) 



(2.18) 



(2.19) 
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and by the second line and ( 2.1 8| ) we derive that 

1 — z kl+1 ~ 



H — a l+l~ 



1 - z 



Thus we get explicitly 

1 — z - LJ - 1 — z 



1 %-f 1 1 - z^+i-^+i 



l,...,N-l. 



1 = 1,.. .,N, 



(2.20) 



(2.21) 



3=1 



and then by ( |2.17| ) 



1, 

1 — ^i+i— "H- 1 
1 - z 



> k 



N-l 

n 

2=i+l 



1 — 2 fe3 + 1_ ' C3 ^ 1 

1 - z ! 



/V 



k l+1 >m>k t 



(2.22) 



where / = 0, . . . , N — 1, and fco = — oo. In fact we got the first line for m > fcjv + 1. 
Value m = k,N was added as from the second line we have that Xk N = 1. In the same 
way the interval of validity of the second equality was extended: we added points m = ki 
as values at m = ki and m = fc;+i differ by the shift I — * I + 1 only. 

Asymptotic behavior at tu-infinity, i.e., at z — > enables us to reconstruct potential, 
as by ( [2.11 ) we have 



j=m-\-l 



On the other side by (|2.22|) 
1, 

l + z(N-l-S mM+i )+o(z), 



z->0. 



m > k N 

k l+1 > m > fc ; 



(2.23) 



(2.24) 



that coincides with (2.23) 



It was already me ntione d above, that only one nontrivial solution exists in t his c ase, 
i.e., solution gi ven in ( 2.22). Indeed, if we construct solution x(z) by means of (2/7) we 
get thanks to (2.8) and ( 2.22 ) relation 



= 9(m > k N )(l - l/z)x 



(2.25) 



I t is also impossible to introduce the Jost solution normalized at — oo, cf. ( |2.9| ), as solution 
( ^.tj ) obeys homogeneous equation 



El — z 



L 0,l 



k-l 



k<l 



X, 



0,k+l- 



(2.26) 



Jost solu tion (2.22) of the equation (2.2) is polynomial with respect to z (i.e., w 



because of (2.5)) with coefficients belonging to N. It is convenient to introduce new 
parametrization of this solution that will be useful in the construction of the integrals of 
motion below. So we introduce discrete measure 



i-l i-1 

M m )= In II ( 1 - 1n+j ) 1n+i + Y[ ( 1 ~ ^m+j ) Qm+i . » 

n>m+l j — 1 j — 1 



(2.27) 
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where m 6 Z and i = 1,2,.... Writing down values for the lowest i's 

fl( m )= E Wn+l+9m+l> ( 2 ' 28 ) 
n>m+l 

/2H = X! ?n( 1_ «n+l)?ri+2 + ( 1 - , Wl) , W2 = ( 2 ' 29 ) 
n>m+l 

./aM = 51 9n( 1 ~9n+l)( 1 -5«+2)9n+3 + ( 1 -9m+l)( 1 -9m+2)9m+3> ( 2 - 30 ) 
n>m+l 

etc., we see that for every m and t first term in f%(m) is equal to the number of higher (i — 
l)-defects, i.e., for all n > m + 1 number of combinations {q n , . . . , q n +i} = {1, 0, . . . , 0, 1} 
with exactly i — 1 zeroes inside, while the second term is equal to 1 iff fc; — m = i for the 
lowest fcj G /C, fc; > m. In particular, we see that if > m > k\ then by ( |2.27 ) 

/<(m) = /i(fc, +1 ) + «fc, +1 - mj i. (2.31) 
Taking into account that by (|1 . 5|) = 1 we get 



= E 9n - W^i, ( 2 -32) 

">fcj i=i 

in particular 

i-l 

M k l) = E Qn Y[( l - Qn+Mn+i = E ( 2 ' 33 ) 
n j=l j>l 

i.e., /i(fci) is the number of neighbor cites in the support K, of g n occupied by units, 

/i(fci) are numbers of (i — l)-dcfccts in the support. As support JC is bounded (see (1.5)) 
we have that 

/ i (fc 1 )=0, i>k N -k 1+ 2, 1 = 1,. ..N, (2.34) 



so that by ( |2.3l| ) 

fi(m) = S kl - m ,i, i>k N -k x + 2. (2.35) 
Also by ( p. 27 ) it is easy to see that 

/■(m)=0, for m>k N , (2.36) 
/i(fc w -l) = l, / 4 (fc Ar -l)=0, *>2 (2.37) 

Definition ( |2.27 ) can be rewritten in the form 

i — 1 i 

/iM = EIl^ 1 9n+j)ln+i E II^ Qn+jhn+i+l (2-38) 
n>m j — 1 n>m j — 1 

that gives the following relations: 

oo 

E^ TO ) = E ?« ( 2 - 39 ) 

i= 1 n>m+l 

y ! /,H = i,' ^S™'- (2 40) 

^ "" v 1 \k N ~m, k N <m, y ' 
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To prove these equalities, notice first that sums are convergent due to (2.35). Then ( |2.39 ) 
follows directly from (2.38), as well as (2.40) for m > k N . Let now k N — 1 > to, then 
thanks to (|2.38|) 



oo i—1 



{k N — l—m i k N -m j 

i—0 j — 1 i— j—1 



that proves (|2.40[) . 

Let u s turn now to ( |2.22| ) for k l+1 > m > fe;, I = 0, . . . , N — 1, ko = — oo. Because 
of (2.33) this equality can be written for k l+1 > to > fc ; as 



n i 



1 - z 



k+l \ a k, 



00 /I /»=( fc !+l) 

n ^ 



fc=i 



Taking ( ^.31 ) now into account we get finally 



n 



fe=i 



(2.41) 



that thanks to ( |2.36 ) gives also valid result in the first line of ( 2.22| ) at m > k 



3 Solutions by mod 2 



It was already mentioned that solution (2.41) is polynomial in z with coefficients belong- 
ing to N, so it solves (2.4) by modulo 2, i.e., it obeys 



x rn-l + (1 + Z)x m + Z(l + q m )x rn+1 = 0, 



mod 2. 



(3.1) 



In [13] it was proved that islands, i.e., intervals separated by at least three zeroes, evolve 
independently. Thus it is natural to consider here the case where support K, is just one 
island. This means that we impose condition that all kj G JC obey 



fcj+i kj < 3, j — 1, 



,N-1, 



or equivalently, 



(! - SnX 1 - - In+i) =0, k N > n > k x 



Then by (2.27) 



/^(m) = 0, i > 4, m>k 1 



(3.2) 



(3.3) 



(3.4) 



so in the case of only one island our solution ( 2.41 ) is given as product of powers of three 
polynomials: 

x m = {1 + zY l(m \l + z + z 2 ) f ^ m \l + z + z 2 + z 3 ) f3( - m \ m>fc 1 -2. (3.5) 

Polynomials 1 + z and 1 + z + z 2 are independent on F2 , but for the third polynomial we 
can write 



1 + z + z 2 + z 3 = (1 + z y 



3d 2, 



(3.6) 
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and because of this 

x m = (1 + z ) h{m)+5f ^ m \l + z + z 2 ) h{m \ mod 2. 



(3.7) 



Moreover, by ( f2~40| ) and Q 

fi{m) + 2/ 2 (to) + 3/ 3 (m) = k N — m, > m > k\ — 2, (3-8) 

so finally we get 

x m = (l + z) fc "- m (l + z)- 2/2 ( m )(l + z + z 2 ) /2 ( m \ mod 2, (3.9) 

where again fc/v>m>/ji— 2. Thus the Jost solution is parametrized by the only one 
dynamical variable /2. 



The original potential can be reconstructed by means of f 2 . Indeed, by ( 2.29 ) 
/ 2 (m - I) - f 2 (m) = (1 - g m )(g m+ i - (1 - q m+1 )q m+2 ) 
that thanks to ( |3 . 3|) can be written in the form 

f 2 (m-l)-f 2 (m) = (l-q m )(2q m+1 -l), k N > m > k t - 2, (3.10) 
and then 

<? ro = l-|/ 2 (m)-/ 2 (m-l)| (3.11) 
as \2q m+ i — 1| = 1. Let us emphasize that equality ( |3.1l| ) is exact, i.e., not the mod 2 



4 Time evolution and integrals of motion 

In order to study time evolution by means of the above construction we need to consider 
two neighbor values of time only, say, t and t + 1. So let us simplify notations used in 
Introduction by denoting 



_ t t 

"n 9n) x n X n1 



etc., 



in 



etc. 



Then equation of motion (1.1) takes the form 

1m = 1m + 9m-29m+l + 9m-l9 m +2' 



mod 2, 



(4.1) 
(4.2) 



(4.3) 



and corresponding transformation of the Jost solution is given due to (1.1) by means of 

mod 2, 



or using notations (2.3) and (2.5) we have 

x m ~ x m ' z< lm-l c lm+2 X m+2i 



mod 2. 



(4.4) 
(4.5) 



Time evolution of /^(tti) is given by substitution of x m in (4.5) by (3.9) 



(1 + z) 2 (W2M)(l + z + z 2 )^ = (1 + z f^S^)){\ + z + z*)Mm) + 

+zq m _ l q m+2 {l + z)- 2 ^ m+2 \l + z + z 2 )^ m+2 \ mod 2. (4.6) 



8 



Considering this (mod 2)-equation separately for the cases q m -iQm+2 equal to and to 
1, we arrive to the exact (not mod 2) equation 

/2M = /2M + 9 m -i9 m + 2 ( 2 9™+i - !), m > fci - 2 (4.7) 

(fc^r > m > ki — 2 is interval of validity of ( |3.9| ) and validity on the interval m > fcjy + 1 
follows from (p6|)). Using ( |3TTo| ) for / 2 we get 



/ 2 (m - 1) - ,/ 2 (m) = 

= (}- 1'n - $ m -lVm+2 + in- 2 9m+l( 2( 7 m ~ !)) ( 2 9m+l ~ *)> 



(4.8) 



and then it is easy to check that I/2 (m — 1) — / 2 (m) | is equal to or 1 , so we can use ( |3.11 ) 
to define 

? m = l-|/ 2 (m-l)-/ 2 (m)|. (4.9) 

Thanks to this relation we der ive the exact eq uatio n of time evolution of q m , that sub- 
stitutes (mod 2)-equation (O). Indeed, using (4J3) we get 



in = 1 - I 1 - 1m - Qm-llm+2 + Qm-2Qm+l ( 2 <?m ~ 1)1- ( 4 -!0) 
The r.h.s. can be simplified if we notice that 

|1 — 1m ~ lm-llm+2 + Qm-2Qm+l ( 2 9m — 1)1 = ^ — 9m + l?m-l3TO+2 — <7m-2 < Zm+l I ( 2 9Vn — 1); 

so finally we get equation 

Qm = 9m — l<?ro-l9m+2 — <7ra-29ro+l I ( 2 9m — !)• (4-11) 

In order to consider integrals of motion we introduce elements of the monodromi 
matrix by means of the standard relation 



x rn =a(z) + z m b(z), m < k 
and analogously for x m . Thanks to (fy*) for m < k\ 



i- 



(4.12) 



thus 



a(z) = a{z), 



(z) = b(z) 



(4.13) 



i.e., bot h a(z ) and b{z) are integrals of motion as it was demonstrated in [13]. Here be- 
cause of ( |3.9| ) we single out functionally independent integrals explicitly. Indeed, by ( [2.17 ) 



b(z) = z 1+kl a(z), (1 - z)a(z) = x kl (z), mod 2, 
so Xfcj (z) is also integral of motion and by ( |3.9[ ) we have 

x kl (z) = (l + z) kN - kl {l + z)- 2f ^ kl \l + z + z 2 ) f ^ kl \ mod 2. 



Now (4.13) and (4.14) are equivalent to 

™i = "'li 



/ 2 (fci)=/ 2 (fci), 



(4.14) 
(4.15) 

(4.16) 
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where the first line was already known (see ( |1.6| )). 

Thus in addition to two known integrals — borders of the island, k% and k^, — we 
proved that /2(fci), number of the single defects (isolated zeroes, see comments af- 
ter ( 2.33| )) inside the island, is preserved under time evolution. This new integral is 
explicitly demonstrated o n Fig . 1. It is cl ear that these three integrals are mutually in- 
dependent and thanks to ( 4.14 ) and ( |4.15| ) they uniquely determine a(z) and b(z). Thus 
all coefficients of the Taylor expansion of a(z) are funct ions of these three integrals only. 
We emphasize that in spite of the fact that Eq. ( 4.1 5| ) is equation by mod 2 all these 
three integrals are exact ones. 

An example of (mod 2)-integral is provided by N. Indeed, in analogy with (3 
have from ( 2.39 ) that 

fx M + f 2 M + /a M = q r 



we 



j>m+l 



On the other side (|3.8|) gives 



/1H+/3H = k 



N 



mod 2, 



so that 



Then 



X] <?j = ^ -™ + / 2 ( m ): mod 2. 

j>m+l 



N= Qm + 1 = k N- k i + 1 + h(ki)i mod 2, 

m>ki +1 

and as any exact integral of motion is also (mod 2)-integral we obtain that 



N = N, 



mod 2. 



(4.17) 



Explicit formulas for /i(fei), /2(&i), and f^(k\) are given in ( 2.33| ) and Fig. 1 shows 
that indeed neither /i(fci), nor fs,(ki) are preserved under time evolution. Of course, 
this does not mean that there cannot be some other integrals of motion. Moreover, it 



is unclear if existence of three exact integrals ( 4.16 ) and one (mod 2)-integral ( 4.17 ) is 
enough for the system to be integrable in some sense. 



One of results of the above construction is substitution of (mod 2)-equation (1.1) 
(or (4.3)) by ( 4.11 ). The latter is the exact one, as for q rn with range equal {0,1} 
we get q m with the same range {0, 1} without using summation by mod 2. Thus the 
problem of integrability of this cellular automaton is reduced to the problem of existence 
of an (exact!) Lax pair for the equation (4.11). This problem, as well as problem of 
definition of the reasonable spectral data for the linear operator (1.1) on F2 are still 
opened. Moreover, we have to stress, that such notions as "integrability" and "degree of 
freedom" themselves must be modified to be applicable to the case of systems on finite 
fields. 
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BBBBBi II BBBB BBl II II BBl BBl II IBBBB 

□□□■■□■□□■■■■□I BBBi II IJBBBIJIJBBDB II Bl II BBI 

Bl BBBi II iBBBBi n MB BBBBi n Bl 1BBBBB1 n 1BB1 n 1BB1 1BB1 n BB 

II II IBBBI IBBI II IBBBBI II IBI II II IBI II IBBBI IBI II IBBBBI II IBBI II IBBI IBBI II IBI II II 
□□□BBDBBBBDDBBBBBB BBl II BDBBBDDBBBBI n BBi 

□□□■□■□□■■Bl B BB BBBBi Bl n BBi II BBBBI II BBI II BBBBLIBrirO 

Bl BBi II BBBBBBi n B Bl n Bl ■■ n ■■ n BBBBi n ■■ n ■■ BB 

BBi BBi II Bl II BBBBB BBBi BBBBi n BBi n BBBBi n ■■ n ■■ B 

Bl BBBBBBBi n Bl n B BBi Bl n BBBi n ■■ n BBBBi n BBi n Bl B 

□□□BBBBBBBLIBBBBBBB Bl BBBIJIJBBBIJIJBBIjriBBBBririBBBriBBriljri 

□□□BLiriBLiriBBBBBBBriB BBBQBBQQBBBQQBBQQBBBBQQBQBBBQQQ 

BBBBBi II Bl II BBi BB BBi BBBBi n BBBi n ■■ n BBBBBi Bl n B 

Bl II BBBBBBi n BBi B Bl Bl n BBBi n BBBi n BBi n Bl n Bl BBBB 

□□□BBI II Bl II BBBi II Bl B Bl IBBIJIJBBBIJIJBBBIJIJBBBBBBriBririBBIjriri 

II II BBBBBBI II IBBBBI IBBI II II IBBI IBBI II IBBBI II IBBBI II IBI II IBI IBBBI II IBI II II 

Bl II BBBBi II BBi BBB Bl BBBBi n BBBi n BBBBBBBi Bl n BBBB 

□□□BBririBBBBririBBBDB BBBBLIBBIjniBBBl n Bl II BBLIBBBririBBBLiriri 

□□□BBBririBBBBririBLiBB bi inBBriBBijriBBBBBBijriBBBriBBririBBririri 

□□□BBBBriLIBBBBBLIBBB BBI II BBI BBl II Bl II BBBI II Bl BBBBI 

II II IBI II IBBI II IBI II IBBBI IBI II II IBBBI II IBBI IBBBBBI II IBBBBI IBI II IBBBBBI II II 

BBl II BBBBBi n Bl BB BBBBi n BBBBBi BBl n BBl BBBi n Bl n B 

BBBi II Bl II BBBi BBB Bl n BBl n Bl n BBl BBl n BBBi BBBBBBB 

□□□BBBBBBBI n ■ ■ B BBl II BBBBBLiriBBriBBririBriB II Bl II BBLILILI 

II II IBI II IBI II IBBBI IBBBBBI II II IBBBI II IBI II IBBI II IBBI IBBBI BBBBBBI II IBI II II 

BBBBBi II Bl Bl n BBB BBBBBBBi n BBl n BBBi BBBBBBi BBBB 

□□□BriLIBBBBLIBBB n BB Bl lDBDDBBI BB 

□□□■■□□■■□■i BB B ■■BBBririBBririBBBriBBBBBBririBBLIBBLiriri 

BBBi II Bl BBl II BBBB Bl n BBBi n BBl n Bl Bl n Bl n BBl n BBl B 

BBBBBi BBBBi n BBB BBl n BBBi n BBBi BBBBBBi n BBl n Bl B 

Bl II Bl Bl II BBBi II BB BBBi n BBBi n Bl Bl n Bl n BBl n BBBi BB 

□□□BBBIJBBBLI BBBI II B BBBBLiriBBBBLIBBBBBBri BBI II 

BBl Bl II BBl II BBBBB Bl n BBl n BBl Bl n Bl n BBl n BBBi Bl n B 

Bl BBBi II BBl II Bl n B BBl n BBl n Bl BBBBBi n BBl n Bl BBBBB 

BBBi BBl II BBBBBBB BBBi n BBBi Bl n BBBBi n BBBi Bl n BBB 

□□□BBriBBBBIjriBljriBBB BBBBI n Bl BBBijriBBBBririBriBBBririBBriro 

Bl BBl II Bl II BBBi n B BBl n Bl BBBi n BBl n Bl n Bl BBl n BBBB 

BBl BBBBBi n BBBBB BBBBi Bl n BBl n BBBBBBi BBBBi n BBB 

Bl Bl II BBBBi II Bl n B Bl n Bl BBl n BBl n Bl n Bl Bl n BBBi n BB 

□□□BLIBBririBBBBBBBBB BBBI BBBBLiriBBBBBBLIBBBI I BBBI 

BBl BBl II Bl II Bl II BB BBl Bl n BBBi n Bl n Bl Bl n BBl n BBBBB 

Bl ■■■■■■■■■■ II B Bl BBBi II BBBBBBBi BBBi n BBl n Bl n B 

BBBBBBBBBBi BBBB BBBi BBl n Bl n BBl Bl n BBl n BBBBBBB 

BBLIBBBBBBBLir ■ BBI II BBI II Bl II 

□□□BBBBBBBBIJIJBBIJBB B B Bl n BBBBI BBBBI n BBBBBIjriBBririri 

BBBBBi II BBBi n Bl B BBBBBi BBl n Bl BBl n BBBBBBi n BBBB 



Fig. 1. Here the horizontal axis is the n-axis and the downward vertical axis is the t-axis. 
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